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A mixed Heisenberg spin chain with frustrated side chains is investigated by numerical and per- 
turbational calculations. A frustration-induced quantum partially polarized ferrimagnetic phase 
and a nonmagnetic spin quadrupolar phase are found adjacent to the conventional Lieb-Mattis type 
ferrimagnetic phase or the nonmagnetic singlet cluster solid phases. The partially polarized ferri- 
magnetic phase has an incommensurate spin structure. Similar structures are commonly found in 
other frustration-induced partially polarized ferrimagnetic phases. Numerical results also suggest a 
series of almost critical nonmagnetic ground states in a highly frustrated regime if the side chain 
spins weakly couple to the main chain. 
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I. INTRODUCTION 

The interplay of frustration and quantum fluctua- 
tion has been extensively studied in a variety of low- 
dimensional quantum magnets. Even in one-dimensional 
cases, various exotic quantum phenomena such as spon- 
taneous dimerization, [1] 1/3-plateau with spontaneous 
trimerization, [2] and transition between quantum and 
classical plateaus [3] are reported. On the other hand, 
the mixed quantum spin chains also have a variety of 
ground states ranging from quantum ferrimagnetism [4] 
to spin gap phases [5-8]. 

Recently, it has been reported that frustration induces 
a partially polarized ferrimagnetic (PPF) phase [9-12] in 
addition to a conventional Lieb-Mattis type ferrimagnetic 
(LMF) phase. The PPF phase appear when both frus- 
tration and quantum fluctuations are fairly strong. It is 
an interesting issue how general such a ferrimagnetism is. 
Hence it is important to investigate the features of PPF 
phases in various spin systems. We are then motivated to 
find a PPF phase in other models and investigate them 
in detail. 

We have introduced a spin chain with side chains in a 
previous paper [8]. This spin chain has frustration ow- 
ing to the interaction among spins in the main chain and 
those of side chains. The frustration varies in strength 
and in feature with the variation of parameters in the 
model. Since we focused on the spin-gap phases in [8] , we 
have investigated the parameter regimes where frustra- 
tion is not strong enough to destroy the spin-gap phases. 
We then found two spin-gap phases and explained them 
by singlet cluster solid (SCS) pictures. 

In the present work, we examine this model in a highly 
frustrated regime. This regime, in which the frustration 
plays a central role, is of interest in its own right, since 
the model exhibits features very different from those in 
the weak frustration regimes. We actually find clear nu- 
merical evidences not only for the above mentioned fas- 
cinating PPF phase but also for the spin quadrupolar 



(QP) phase [13-19]. These phases are totally different 
from the conventional phases such as spin gap phases and 
the LMF phase which can be realized even in the unfrus- 
trated case. They will be investigated in detail in the 
present paper. Also numerical data are obtained which 
suggest the possible existence of an exotic almost critical 
nonmagnetic ground state in the regime where the cou- 
plings between the side chain and main chain spins are 
weak but strongly frustrated. 

The transition from a ferrimagnetic phase to a non- 
magnetic phase in the present model takes place because 
the quantum fluctuation in the side-chains destroys the 
ferrimagnetic long range order in the main chain. The 
mechanism of quantum destruction of ferromagnetism 
and ferrimagnetism has been less studied than that of 
antiferromagnetism; the latter has been extensively stud- 
ied in relation with the high-Tc superconductivity. Re- 
cently, however, experiments have been reported on the 
nonmagnetic ground states in one- and two-dimensional 
materials [23, 24] with ferromagnetic nearest neighbour 
and antiferromagnetic next nearest neighbour couplings. 
Theoretical investigation has also been carried out for 
corresponding models [25-28]. 

In an unfrustrated ferrimagnet, the spontaneous mag- 
netization is uniquely determined by the Lieb-Mattis the- 
orem [29]. This type of quantum ferrimagnetism has 
been investigated in detail [4]. As far as the frustration 
is weak, the spontaneous magnetization remains locked 
to this value [30]. This phase is the LMF phase [10]. 
The spontaneous magnetization in this phase is a sim- 
ple fraction of the saturated magnetization. In contrast, 
the spontaneous magnetization in the PPF phase con- 
tinuously varies with the parameter characterizing the 
strength of frustration and is not a simple fraction of 
the saturated magnetization. The PPF phase appears 
between the LMF phase and the nonmagnetic spin gap 
phases. This type of phase is first predicted in the pio- 
neering work of Sachev and Senthil [21] in the quantum 
rotor model. Bartosch, KoUar and Kopietz [22] proposed 
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a possibility of ferromagnetic Luttingcr liquid in an itin- 
erant one-dimensional Fermi system. The first explicit 
example of quantum PPF phase induced by frustration 
in one-dimensional quantum spin systems was proposed 
by Ivanov and Richter [9] in a frustrated mixed spin lad- 
der. Similar phases are also found by Yoshikawa and 
Miyashita [10] in a uniform spin chain and by one of the 
present authors in a trimerized zigzag chain [11, 12]. We 
propose another example of the PPF phase in the present 
mixed spin chain. The present example is substantially 
different from previous ones, because it is accompanied 
by the destruction of the ferrimagnetic order in the main 
chain by the frustrated coupling to the quantum fluctu- 
ation in the side chains. 



The QP phase is well known for the spin-1 bilinear- 
biquadratic chain between the Haldane and ferromag- 
netic phases [13-19]. The exact Bethe ansatz solution 
is available if the coefficients of the bilinear and the bi- 
quadratic terms coincide with each other [13-15]. Re- 
cently, similar phases are found in the frustrated two- 
dimensional Heisenberg model with ferromagnetic near- 
est neighbour interaction and antiferromagnetic next 
nearest neighbour interaction [28]. In this paper we 
explicitly show that our model reduces to the bilinear- 
biquadratic chain in appropriate limiting cases. It is also 
argued that the QP phase should appear in a wide class 
of complex spin models between ferrimagnetic and spin 
gap phases. 



This paper is organized as follows. In the next section, 
the model Hamiltonian is presented. Various limiting 
cases are discussed using the perturbational approxima- 
tion from the strong coupling limit in Sec. III. The nu- 
merically obtained ground state phases are explained in 
Sec. IV. The properties of PPF phase are described in 
detail in Sec. V. Section VI is devoted to summary and 
discussion. 



II. HAMILTONIAN 

We consider the mixed Heisenberg spin chain described 
by the Hamiltonian 

iV/3 

'H = Y, [JlSi(p)S2(p) + J2S2(p)Si(p+ 1) 

+i^iSi(p)T(p) + K2S2(p)T(p)] . (1) 

where ^2{p) and T(p) are the spin-1/2 operators and 
Si(p) is the spin-1 ones in the p-th unit cell, as shown 
in Fig. 1. In what follows, we use nondimensional pa- 
rameters j = J2/J1, k = K\l J\ and r = K^ji^Kx), and 
the unit of Ji = 1. The parameter r characterizes the 
strength of frustration. The total number of spin sites is 
denoted by TV, and then the number of unit cells is iV/3. 

In regime < r < 0.5, where frustration is not strong, 
we have found two types of nonmagnetic ground states 
and have explained them by SCS pictures [8]. When 
frustration becomes strong, the phase diagram drastically 
changes. In the present paper, we will investigate the 
strongly frustrated case r > 1 in detail. 



III. LIMITING CASES 

A. Small j regime 

For i = 0, the chain is decoupled into an assembly of 
3-spin clusters described by the Hamiltonian 



n^ip) = Si(p)S2(p) + fcSi(p)T(p) + 2krS2ip)Tip). 



(2) 



As discussed in Ref. [8], the cluster ground state is a 
singlet state with energy 
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for k < kr 



1, and a triplet state with energy 




Si(p) S2(p) 



FIG. 1: The quantum spin chain with side chains which we 
study in this paper. Si(p), S2(p) and T(p) are spins in the 
p-th unit cell. These magnitudes are S\ = 1, S2 = \, and 



T = i, respectively. 



-(1 + fc + 2kr) - y/(l + k- AkrY + S{k-iy 



(4) 



for k > kc- Therefore, the ground state of the chain for 
small j is a gapped local 3-spin singlet phase for k < 
fee. We call this phase as Gap I phase following Ref. [8]. 
Because this cluster ground state is gapped, it cannot 
gain energy within the first order in j even in the presence 
of j. 

In the 3-spin triplet ground state for fe > fee, we define 
a composed spin with magnitude 1 as S{p) = Si{p) + 
^•lip) + T(p) to describe low energy phenomena. The 
total Hamiltonian is written as Ti. 



Hf^^ with an 
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unperturbed part Ti.^ 
part 



^ H.'^ {p) and an interaction 



QP 



njA 



N/3 

= 5^iS2(p)Sib+l). 

p=l 



(5) 



The perturbation calculation up to the second order in 
the interaction part H^-^^ yields the following effective 
bilinear-biquadratic Hamiltonian for the composed spin 
Sip): 



njA 



N/3 

J2 [je1tS(p)S(p+ 1) + D^^{S{p)S{p+l))' 



(6) 



where the effective interaction parameters consist of the 
first and the second order perturbation terms as — 



D 



A(l) 



D. 



A{2) 



In the first order perturbation, the effective interaction 



^(1) 



parameters coming from J'^^ 
D. 



and D^^^ are 



jA 
off 

-.A 

off 



jX{a), 



0, 



(7) 

(8) 



where X{a) is given by 
X{a) 



with 



\/2(l + a2)2 V 2^/2 



2V2(1 - fc) 



^{k + 1 - Akrf + 8{k - 1)2 -l-k + ikr 



(10) 



For small but finite j, the energy of the LMF phase is 
given by Et + per unit cell while the energy of the 
3-spin singlet (Gap I) phase is given by Es with no first 
order correction in j. Therefore, comparing the energies 
of these two ground states, we find that the phase tran- 
sition between the gapped 3-spin singlet phase and the 
LMF phase takes place at 



k = /stf =kc 



(l-r)(2-2r + r2)(3-r)j 



2r(2-r)(3 

7^(1) 



4r + 2r2) 



(11) 



The effective coupling J^^^ vanishes for k = 1. There- 
fore, in the neighbourhood of fc ~ 1, the higher order 
terms come into play. Within the second order pertur- 
bation with respect to 1 — A; and j, the effective coupling 
constants are 



J. 



J 



eS 



J-(l - fc) 

8(r-l)(l + 2r) 2(2r - 1) ' 



8(r- l)(4r2 - 1)' 



(12) 
(13) 



The effective model (6) has a variety of phases [19]. 
Within the present parameter regime, we find the Hal- 
danc phase for < < which corresponds to Gap 
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FIG. 2: Phase diagram for small j with r = 1.2. 

II phase in Ref. [8], the QP phase for < < 
and the LMF phase for < 0. However, in the orig- 
inal Hamiltonian (1), the LMF phase is limited by the 
transition to the 3-spin singlet phase at = Es — Et 
as discussed above. For > Eg — Et, the ground state 
is the 3-spin singlet (Gap I) phase. 

Thus the conditions for each phase in terms of the orig- 
inal parameters are summarized as follows: The ground- 
state phase of the present model is (i) the 3-spin singlet 
(Gap I) phase for 



< < fcxF, 
(ii) the LMF phase for 



/ctf < k < fcpQ — 1 — 



(2r - l)j 



4(r-l)(2r+l)' 



(iii) the QP phase for 



fcpQ < k < fcoH = 1 — 



QH 



J 



2(2r + l)' 

and (iv) the Haldane (Gap II) phase for 
k > /cQH. 



(14) 



(15) 



(16) 



(17) 



These phase boundaries are plotted on the k-j plane in 
Fig. 2 for r = 1.2. It should be remarked that the 
QP phase in the present model is realized without bi- 
quadratic interaction in the original Hamiltonian (1). 



B. Large j Regime 

In the large j limit, Si(p) and S2(p— 1) form an effec- 
tive S = 1/2 spin (t{p) = Si(p) + S2(p - 1). The effective 
Hamiltonian is given by 

N/2 

H = Y,{Kfa{p)T{p) - Jfa{p)a{p + 1) 



-Kfcr{p+l)T{p)), 



(18) 
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FIG. 3: A-chain realized in the limit j >> 1. 



which form a A-chain structure depicted in Fig. 3. The 
effective interactions are given by 



9' ^ 



3 ' ^1 



as argued in Ref. [8] . The detailed analysis of this model 
is reported in a separate paper [20]. Therefore, we only 
quote the results and rewrite them in terms of the original 
model (1). 

The ferromagnetic phase of the model (18) corresponds 
to the LMF phase in the original model (1). This phase 
is stable for 



< fc < fcpQ = 



r-2 
3r 



(20) 



even in the limit of large j. 

For k > fcpQ, the ground state is nonmagnetic. Never- 
theless, there are still several different phases. For large r, 
the model (18) reduces to the S ~ 1 bilinear-biquadratic 
model and the QP phase appear for 



fcpQ < fc < fcoH = 



3r 



QH 



9r 



(21) 



and Haldane phase appear for fc > fcqH- 

The results of the numerical diagonalization calcula- 
tion for the model (18) is summarized in Fig. 4. There 
are the QP, the Haldane and the LMF phases, which we 
discussed above. In addition, numerical results suggests 
that there possibly exist a narrow PPF phase between the 
LMF phase and the QP phase, and almost critical non- 
magnetic ground states for small values of fc. The PPF 
phase is so narrow that it cannot be represented in Fig. 
4. We speculate that the almost critical nonmagnetic 
phases are spin gap phases with extremely small energy 
gap with large scale resonating singlet cluster solid struc- 
ture. Corresponding phases are also found for finite j as 
described in the next section. 



C. Large r Regime 

We examine the case of K2 ^ i^i, Ji (r, fcr ^ 1) in 
this subsection. In the the limit oi r,kr —^ 00, the chain 
is decoupled into an assembly of 3-spin clusters, each of 
which described by the Hamiltonian 

H''ip)=jSi{p+l)S2{p) + 2krS2{p)T{p). (22) 
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FIG. 4: Phase diagram in large j limit. The shaded region is 
the almost critical nonmagnetic phase. 



The eigenvalues of this 3-spin Hamiltonian are given by 



E{S = l;g) 
E{S=l;e) 



kr j 

—j — 2kr — y^l6k'^r^ — 8krj + 



-j - 2kr + v'lefc^r^ - 8fcrj + 9f 



E{S = 0)^--j. 



(23) 
(24) 

2 

(25) 
(26) 



The ground states are the triplet states with energy 
E{S = 1; g). It should be noted that the lowest excitation 
energy is of the order of kr even if j is small. Therefore 
the perturbation calculation from this limit is valid even 
for small j. As a result, each cluster has an effective spin- 
1 degree of freedom S{p) ee Si(p+1) + S2(p) + T(p). The 
total Hamiltonian is written as 7i = "Hq + Ti-^i with an 
unperturbed part TLq = Ti^ [p) and an interaction 
part 



N/1, 

W?nt=E(S2(p) + fc-T(p))Si(p). 

p=l 



(27) 



We can write down the effective Hamiltonian for S(p) up 
to the second order in Ti.^^^ in the form. 



Af/3 



^ [ji|S(p)S(p + 1) + 2?fff (S(p)S(p + l)f 



(28) 



Within the first order in Ti^j, the effective interaction 
parameters are given as 



eff 



-{X{a) + kX{-a)), 



D. 



eff 



0, 



(29) 
(30) 
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FIG. 5: Phase diagram on the j/r-rAfc plane where Afc = 
k - kc{j/r). 



where X{a) is given by Eq. (9) with 

7 - 4kr + Jiek^r^ - 8krj + 9f 

a = — ■ 

2V2j 



(31) 



The effective exchange constant J^g vanishes up to the 
first order in if j and k satisfy the relation 



8k{P - 1) 
(3A;-l)(3fc-5)' 



(32) 



We denote the value of k which satisfies this relation by 
kc{j/r) as a function of j/r. The terms of 0{r^^) come 
into play for k ~ fee (j/r), and constitute the bilinear- 
biquadratic form (28) with ~ 0{k — kc,r~^) and 

We do not explicitly present the second order expres- 
sion for and -D^, since we numerically carried out 
the summation over the intermediate states in the sec- 
ond order perturbation calculation. The LMF-QP phase 
boundary fcpQ is determined by setting = 0. Be- 
cause the correction terms are of 0(r~^), deviations 
A/cFQ = fcpQ - kc{j/r) and A/cqh = fegn - fcc(j/?') 
scale with 1/r for fixed j/r. Figure 5 shows the j/r- 
dependence of rAfepq and rAfequ. The phase boundaries 
for r = 10 determined by the present approximation are 
shown in Fig. 6. 

The calculation in this subsection suggests that the QP 
phase found in the small-j limit (III A) and that in the 
large-j limit (HIB) form a single phase, although it is 
explicitly demonstrated only in the large r limit. 



IV. NUMERICAL GROUND STATE PHASE 
DIAGRAM 

For small r, there are two types of nonmagnetic ground 
states and the Gaussian transition occurs between them 
as described in Ref. [8]. The pcrturbational approaches 
in III, however, predict the presence of the LMF and 
the QP phases in addition to the conventional spin gap 
phases. 

We start with the case of r = 2, where the frustration 
is fairly strong. The ground state phase diagram is shown 



40 



20 





\ r=10 




i-li— QP 




W Gap II 


LMF 


\\ (Haldane) 



0.5 



k1 



FIG. 6: Phase diagram for r = 10 using the approximation in 
subsection IIIC. 
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FIG. 7: Ground state phase diagram of (1) for r = 2.0. In 
this and following figures, phase boundaries determined from 
the numerical diagonalization data for N = 18 and 24 are 
shown. The solid lines are guide for eye. The broken lines are 
the results of the small-j approximation in III A. 



in Fig. 7. The phase boundaries calculated by using the 
numerical diagonalization data for = 18 and 24 are 
shown. Between the Haldane-like Gap II phase and the 
LMF phase with spontaneous magnetization Mq ~ Ms/2, 
there appear a QP phase, which is identified by the low- 
est excitation with total spin 2, [16, 17] for 0.66 ^ A: < 1 
as expected from the pcrturbational calculation. For 
k > 0.5, only the data for iV = 18 are shown consid- 
ering the quasi-trimcrized nature of the QP phase. In 
spite of the limited system size, the phase boundary co- 
incides well with the pcrturbational results for small j 
as depicted by the broken lines. However, the QP phase 
vanishes when k decreases and j increases. Instead, there 
appears a PPF phase with the spontaneous magnetiza- 
tion of intermediate values between Ms/2 and 0. The 
detailed properties of this phase is discussed separately 
in the next section. 
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FIG. 8: Ground state phase diagram of (1) for r = 1.5. 




FIG. 9: Ground state phase diagram of (1) for r = 1.3.. 



For r < 2, the 3-spin singlet Gap I phase appears for 
Q < k < kc- As shown in the phase diagrams Fig. 8 for 
r = 1.5, Fig. 9 for r = 1.3 and Fig. 10 for r ^ 1.2, 
both the QP phase and the PPF phase shrink to regions 
around fc ~ 1 with the decrease of r. At r = 1, the 3-spin 
singlet phase extends up to fe = 1 for small j, and both 
the QP and the PPF phases vanish. 

The phase boundary between the Gap I and the 
Gap II phases is determined by the twisted boundary 
method [31, 32]. In Fig. 8 for r = 1.5, the Gap I phase 
consists of two separate regions which have the same 
parity under the twisted boundary condition. However, 
these two regions merge with the decrease of r as shown 
in Fig. 10 for r = 1.2. Therefore wc conclude these two 
regions belong to a single phase. 
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FIG. 10: Ground state phase diagram of (1) for r = 1.2. 



For small fc, we find the region in which the singlet- 
triplet excitation gap Ai? behaves almost critically. Fig- 
ure 1 1 (a) shows the j-dependence of the scaled gap NAE 
for r = 1.3 and k = 0.12 for N = 12, 18 and 24. Al- 
though the boundary of such region cannot be precisely 
determined, they roughly correspond to the shaded re- 
gions in Figs. 8, 9 and 10. Applying the twist boundary 
method, [31, 32] we calculate the ground state energies 
i?+ and E- with spin inversion parities -I- and — , respec- 
tively. Then we find that the spin inversion parity of the 
ground state, or the sign of E+ — E-, changes several 
times with the variation of j as shown in Fig. 1 1 (b) . In 
this region, the difference E+ — E^ is extremely small 
(typically less than 0(10^"^) for N = 24), and becomes 
even smaller with the decrease of k. This behavior is 
most prominent for j ^ 1 as shown in Figs. 9 and 10, 
unless this regime is covered by the ferrimagnetic phase 
as in Figs. 7 and 8. 

The critical values of j at which the parity changes 
are shown by the dotted lines in Figs. 8, 9 and 10 for 
each system size. They depend sensitively on the system 
size. Due to the limitation of the system size, we cannot 
conclude whether these lines correspond to some phase 
transitions in the thermodynamic limit. However, for 
large j, these lines are expected to be continuously con- 
nected to the similar lines of the effective model (18) in 
the corresponding regime (the shaded region in Fig. 4). 
The numerically estimated values of the central charge 
c of the effective model (18) on these lines suggest that 
they are Gaussian transition lines with c ~ 1 among spin 
gap phases with extremely small gap and large scale sin- 
glet clusters [20]. Therefore it is likely that these lines 
in the present model are also similar Gaussian transition 
lines. However, considering the large ambiguity in the 
estimation of c in Ref. [20] , other possibilities cannot be 
ruled out. The elucidation of the nature of the ground 
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FIG. 11: (a) The j'-dependence of the scaled singlet-triplet en- 
ergy gap NAE with periodic boundary condition and (b) the 
energy difference between the different parity ground states 
E+ — E- with twisted boundary condition. The parameters 
are r = 1.3, k = 0.12 and N = 12, 18 and 24. The triangles 
indicate the values of j where the ground state parity changes 
under the twisted boundary condition for N = 24. 
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FIG. 12: The j-dependence of the spontaneous magnetization 
for k = 0.5 and r = 1.3. The results are calculated in the 
N = 72 system by the DMRG method. The magnetization is 
normalized by the saturated magnetization Ms = 2N/3. 
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FIG. 13: Magnetization curve in the PPF phase with k = 
0.4, j = 1.7 and r = 1.5 for iV = 192 calculated by the DMRG 
method. 



state in this regime is left for future studies. 



V. PARTIALLY POLARIZED 
FERRIMAGNETIC PHASE 

A. Numerical Results 

To clarify properties of the PPF phase, we calculated 
the spontaneous magnetization Mq of the ground state 
by the density matrix renormalization group (DMRG) 
method. In Fig. 12, we show Mq as a function of j 
for r ~ 1.3 and k = 0.5 with N = 72. For small j, 
the ground state is in the Gap I phase with AIq ~ 0. 
When j increases, the LMF phase sets in where Mq = 
Ms/2. The slight deviation of Mq/Ms from 0.5 in this 
phase is due to the boundary effect inevitable for the open 
boundary DMRG. With further increase of j, the ground 
state enters into the PPF phase where the spontaneous 
magnetization gradually decreases down to zero. 

Typical magnetization curve calculated by the DMRG 



is presented in Fig. 13 for r = 1.5, j = 1.7 and k = 0.4. 
The magnetization increases continuously from the zero 
field value in the PPF phase in contrast to the LMF 
phase where the magnetization is quantized to the zero 
field value up to a finite critical field [4]. This implies 
that the magnetic excitation is gapless in the PPF phase 
as in the previously reported systems [9, 21, 22]. 

The local magnetization profile {S^{p)) calculated by 
the DMRG is plotted against p in Fig. 14 for r = 1.5, j = 
1.7 and k ~ 0.4. In addition to the period 2 oscillation, 
an incommensurate modulation is clearly observed in the 
local magnetization. 

Similar behaviors are found in other frustration in- 
duced PPF phases in the spin-1/2 period-3 chain with 
next nearest neighbor interaction [11] and in the model 
of Rcf. [10]. We expect these features are common as- 
pects of the frustration induced quantum PPF phase. 
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FIG. 14: Local magnetization profile in PPF phase with k — 
0.4, j = 1.7 and r = 1.5 for N = 192 calculated by the DMRG 
method. 




FIG. 15: Classical planar spin configuration in a triangle. 



B. Classical Picture 

To understand the physical picture of the PPF phase, 
we consider the classical limit of the Hamiltonian (1). 
Since the J2-bonds are not frustrated, the relative angles 
between the spins Si(p), S2{p) and T(p), which form a 
triangle, are not affected by j in the absence of magnetic 
field. The ground state of the whole chain can be con- 
structed by arranging the triangles so that Si{p) and 
S2{p— 1) are antiparallel. The classical ground state en- 
ergy Eq of a 3-spin cluster consisting of the spins Si(p), 




Si(1) 82(1 )Si (2) 82(2) Si(3) 82(3) 
(b) 

>T(1) iT(2) ^(3) 

i/ov yo\. yo\ 

81 (1) 82(1)81(2)82(2) 81 (3) 82(3) 

FIG. 17: Examples of classical ground state configurations (a) 
with finite magnetization and (b) with vanishing magnetiza- 
tion. The inner arrows indicate the direction of the rotation 
of spins along each triangle. 



M/M. 




k=0.2 

k=1 .0 — 
k=2.0 — - 



H 
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FIG. 18: Classical magnetization curves for k = 0.2, 1.0 and 
2.0 with j = 1.7 and r = 1.5. 



S2(p) and T(p) is given by 



Eq {4), 6*, X) = 2 - ^) + Y cos(^ -X) + ^ cos(x - 9), 

(33) 




FIG. 16: Classical phase diagram on the k-r plane. 



assuming a planar configuration. We confine ourselves 
to the planar configuration, because nonplanar configu- 
rations have higher energy. The angles 9, (j) and x denote 
the polar angle of Si(p), T(p) and S2(p) measured from 
2-axis as shown in Fig. 15, respectively. By minimiz- 
ing Eq , we find that a stable noncoUinear configuration 
within the triangle is realized in the region 



< k < 



(34) 



which is indicated as PPF in Fig. 16. Outside this re- 
gion, three different coUinear configurations are flavored 
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depending on the values of k and r as 



\s^{p)T^{p)s^ip): 



k > 



r + 1 



mi) o<k< 



1 

1 — ?■ 



(35) 



If these are regularly arranged along the chain keeping 
the spins Si (p) and S2{p— 1) antiparallel, the LMF phase 
and two kinds of Neel phases are realized as indicated in 
Fig. 16. 

In the region (34), however, the energy is invariant un- 
der the simultaneous rotation of two of the spins (say 
82(73) and T(p)) around the remaining one (say Si(p)). 
Therefore the ground state of the whole chain has a 
macroscopic degeneracy. Among these highly degenerate 
ground states, the states with various values of magneti- 
zation are included. For example, the state depicted in 
Fig. 17(a) has finite magnetization, while the other one 
depicted in Fig. 17(b) has a spiral structure and has no 
net magnetization. 

In the presence of magnetic field, the degeneracy is 
lifted and one of the ground states which has the largest 
magnetization is selected by an infinitesimal magnetic 
field. In this state, the spin configuration should be re- 
stricted in a single plane, because the tilt out of the plane 
reduces the net magnetization. Then each triangle in the 
spin chain takes one of twofold degenerate spin configura- 
tions with clockwise and counterclockwise spin rotations. 
These two configurations are distinguished by chirality. 
In a uniform array of triangles with the same chirality, the 
net magnetization tends to be averaged out. Therefore, 
the magnetization per triangle decreases with increasing 
length of the array. To maximize the magnetization, an 
array of the triangles with alternating chirality depicted 
in Fig. 17(a) is most favorable. This also explains the 
period 2 oscillation observed in Fig. 14. With the in- 
crease of the magnetic field, the spin orientation gradu- 
ally changes to the direction of the magnetic field. 

The ground state in a finite magnetic field H is ob- 
tained by numerically minimizing the classical energy 



JV/3 p 



^ cos{(t){p) - e{p)) + y cos(0(p) - xip)) 

(36) 

1 j 
+ 2 cos(x(p) - 9{p)) + - cos(x(p - 1) - e{p)) 

(37) 

-i?(cos6'(p) + icos0(p) -f icosx(p)) . (38) 

where the magnetic field is in z-direction. Consider- 
ing the argument in the preceding paragraph, we as- 
sume a period 2 structure in the minimization procedure. 
Typical magnetization curves are shown in Fig. 18 for 
r = 1.5 and j = 1.7 with k = 0.2, 1.0 and 2.0. At first 



sight, these three curves correspond to the LMF phase 
[k = 0.2), the PPF phase [k = 1.0) and a nonmagnetic 
phase (fc — 2.0). However, the classical phase including 
the point of A: = 2.0 correspond to a Neel ordered state 
I 5f(l)r-(l)5|(l)5f(2)T^(2)5|(2)...) = | Wi • • • ) 
as discussed above. Nevertheless, we may regard this 
Neel phase as a classical counterpart of the Gap II (Hal- 
dane) phase, since the short range antiferromagnetic cor- 
relation between the total spins of the 3-spin clusters is 
common in the Neel state and in the Haldane-like state, 
and the latter is more stable in the quantum case. Simi- 
larly, in the classical PPF ground state has noncoUinear 
spin structure with broken U{1) symmetry around the 
z-axis, while in the quantum case, the U{1) symmetry is 
restored due to quantum fiuctuation. It should be also 
noted that another ordered ground state with structure 
I iTiitrii •••) corresponding to the classical counterpart 
of the 3-spin singlet phase appears for small r as shown 
in Fig. 16. 

Thus, all the phases in the quantum model (1) have 
the classical counterparts except for the QP phase which 
is of essentially quantum origin. Therefore the classical 
picture appears to be satisfactory at least qualitatively. 
The parameter regime for each phase is, however, largely 
different from that of the original quantum model (1). 
For example, the condition for the classical noncoUinear 
spin configuration (34) is independent of j and does not 
cover the numerically obtained region of the PPF phase. 
In addition, the incommensurate modulation of the spin 
profile is not explained in the classical model. These 
features are essentially quantum effect which is beyond 
the classical interpretation. 



VI. SUMMARY AND DISCUSSION 

A mixed spin chain with frustrated side chains is in- 
vestigated, when the frustration is strong. Not only the 
LMF phase but also the PPF phase appears between 
the nonmagnetic phase and the LMF phase. The PPF 
phase has continuously varying spontaneous magnetiza- 
tion, which is not a simple fraction of the saturated mag- 
netization. The local magnetization profile has an incom- 
mensurate structure in the PPF phase. These features 
are common with other examples of PPF phases induced 
by frustration. Classical interpretation of the PPF phase 
is also presented. It is pointed out that an infinitesimal 
magnetic field selects the PPF state in an appropriate 
parameter region. For the quantum model, however, the 
PPF state is realized in the absence of magnetic field. 
This suggests that quantum fiuctuation selects one of the 
classical ground states. 

The presence of QP phase is demonstrated by the per- 
turbational calculation as well as the numerical method. 
Although the present perturbational calculation is car- 
ried out for our specific model (1), the derivation of the 
effective spin-1 bilinear-biquadratic chain is quite gen- 
eral and is applicable to models which have an effective 
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spin-1 degree of freedom in each unit cell. The QP phase 
appears around the point where the first order effective 
coupling vanishes due to frustration. Hence the QP phase 
is expected to be commonly found in a wide variety of 
frustrated quantum spin chains. 

For r > 1 and j ~ 1 and A; ^ 1, the main chain couples 
only weakly with side spins but the frustration is strong. 
In this case, there appears an almost critical nonmagnetic 
ground state in which the spin inversion parity under the 
twisted boundary condition changes many times with the 
variation of parameters. Considering the continuity to 
the similar ground state in the effective model for large j, 
it is likely that a series of Gaussian transitions take place 
among gapped phases with an extremely small gap. 

Although the nature of this almost critical ground state 
remains unresolved, we may speculate its physical origin 
in the following way: In the absence of side spins T(p), 
the ground state of the main chain is ferrimagnetically 
ordered. For small k, side spins are coupled to this fer- 
rimagnetic moment antiferromagnetically via Ki bond 
and ferromagnetically via K2 bonds. For r ~ 1, the ef- 
fective coupling is even weakened due to frustration, and 
for J = 1, the ferrimagnetic state of the main chain has 
no local valence bond structure. Therefore the effective 
coupling among the side spins, which is mediated by the 
fluctuation in the main chain, would be very long ranged 
for j ~ 1. This implies that the resultant nonmagnetic 



state should have a highly nonlocal character. Thus we 
may speculate that the ground state has an extremely 
small energy gap and large scale singlet clusters. A sim- 
ilar ground state with extremely small gap is known in a 
S ~ 1/2 zigzag chain with ferromagnetic nearest neigh- 
bor interaction and antiferromagnetic next nearest neigh- 
bor interaction [27]. However we do not flnd an explicit 
mapping of the present model onto the field theory of 
Rcf. [27]. 
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